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In their recent paper [Phys. Rev. A 71, 033622], B. T. Seaman et al. studied Bloch states 
of the condensate wave function in a Kronig-Penney potential and calculated the band structure. 
They argued that the effective mass is always positive when a swallow-tail energy loop is present 
in the band structure. In this comment, we reexamine their argument by actually calculating the 
effective mass. It is found that there exists a region where the effective mass is negative even when 
a swallow-tail is present. Based on this fact, we discuss the interpretation of swallow-tails in terms 
of superfluidity. 
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In their recent work on the band structure of Bose- 
Einstein condensates in a Kronig-Penney potential 
Seaman et al. have reported that when a swallow-tail 
energy loop is present in the band structure, the effective 
mass, which is known to be closely related to the dynam- 
ical instability of condensates in periodic potentials 0] , is 
always positive and the condensate is dynamically stable 
in the lower potion of the swallow-tail. In this comment, 
however, we recalculate the band structure by solving 
the time-independent Gross-Pitaevskii equation and find 
a region where the effective mass is negative in the lower 
portion of the swallow-tail. 

We start from the time-independent Gross-Pitaevskii 
equation, 



2m dx 2 



V(x) + g\V (x)\' 



with a Kronig-Penney potential, 

V(x) = V J2Hx-ja), 

3 



*o(a:) = 0, (1) 
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where a is the lattice constant, Vq is the potential 
strength, /i is the chemical potential, ^o(x) is the con- 
densate wave function, and g is the coupling constant of 
the interatomic interaction. Note that we set h = 1. The 
normalization condition is given by 



\* {x)\ 2 dx = N c 



(3) 



where Nq is the number of condensate atoms per site. 
The energy of the condensate per site is given by 



f(j+l)o 

E = I dx^n 
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The condensate wave function can be written as 
^o{x) = ^/noA(x)e lS ( x \ where A{x) and S{x) mean 
the amplitude and phase of the condensate. A(x) is 
normalized by the density at the center of each site 
n = |* ((j + 1/2) a) | 2 . Thus, Eq. Q is reduced to 



1 d 2 A Q 2 



2m dx 2 2m 



A~ 3 + V(x)A- liA + gn A 3 = Q, (5) 



A2 f = Q- 

dx 



(6) 



Equation ([6]) is the equation of continuity and Q describes 
the superfluid momentum. 

Assuming that the condensate sits in the first Bloch 
band, one obtains the solution of Eq. © in the region 
(j - l/2)a <x<(j + l/2)o, 

A(x)* = (l-0-) 



where 



gn 2 2 V \gn 




FIG. 1: First Bloch band of the energy of the condensate 
E(K) for Vb = lgn a v£av and a = 5£ a v, where rw = Nc/a 
and £av = (mgn^ 
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FIG. 2: Group velocity v s (K) for Vo = l<?n av £ av and a = 5£ a 
where c s = (gn av /m) 1/2 . 




FIG. 3: Effective mass m*(K) for Vb = l<?n av £ av and a 
5^ a v - 



sn(w, Z) denotes the Jacobi elliptic sine function. The 



-1/2 



healing length £o is expressed as £o = ( m 9 n o) 
and xq are determined as functions of Vb, a and Q by 
Eq. ([3]) and the boundary condition at x = ja 
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ja+O 



dA 2 
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4mV A 2 (ja). 
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The phase S(x) can be calculated as 

S(x) - S(ja) = J dx-^. 



(9) 



(10) 



Detailed calculations of the condensate wave function 
Eqs. |0) and (JTUJ) have been shown in Ref. [|[. 

Imposing the Bloch theorem ^q(x + a) = ^o(x)e lKa , 
the superfluid momentum Q can be related to the quasi- 
momentum K as 



Ka 



S 



(-!) 



dx ^ 



(11) 



Substituting Eqs. (J?) and (fTTj) into Eq. Q, one obtains 
the first Bloch band of the energy of the condensate E(K) 
as shown in Fig.Q] A swallow-tail energy loop is present 
at the edge of the first Brillouin zone. 

The group velocity v g (K) and the effective mass 
m*(K) are given by 



v s (K) 
m*{K) 



dE 
W 

\dK 2 ) 



(12) 
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v g (K) and m*(K) are shown in Figs. [2] and [3l respec- 
tively. Reflecting the presence of the swallow-tail, v s (K) 
has a reentrant structure. In the lower portion of the 
swallow-tail, v g (K) increases as K increases until Km, 
which gives the maximum value of the group velocity. As 
K increases further from Km to the swallow-tail edge Kb, 
v g (K) decreases. In the upper portion of the swallow-tail, 
v g (K) increases as K increases from ^ to Ke- It is ob- 
vious in Fig. [3] that the effective mass is negative in the 
region of the lower portion between Km and Ke- This 
fact suggests the dynamical instability of the condensate 
in the region, which has been shown in Ref. [3j by the lin- 
ear stability analysis when the lattice constant is much 
larger than the healing length. The region where the ef- 
fective mass is negative is present also in the case of a 
sinusoidal periodic potential, because v g (K) exhibits the 
same reentrant behavior It is worth mentioning that 
the region of the negative effective mass exists for any 
values of Vb and a except in the limit of Vb — > 0. 

In Refs. [3,0], the presence of swallow-tails in the band 
structure was interpreted as a manifestation of superflu- 
idity of condensates. This interpretation is inconsistent 
with the negative effective mass in the lower portion of 
the swallow-tail for the following reason. It is assumed 
in the interpretation that the Bloch states in the lower 
portion correspond to local minima in the energy land- 
scape. However, they do not actually correspond to local 
minima, but to saddle points when the effective mass is 
negative. Moreover, it was insisted in Refs. [l], Q that as 
a superfluid, the condensate can screen out the periodic 
potential and the system follows the energy spectrum in 
the absence of the periodic potential until the superfluid 
momentum reaches the swallow-tail edge. This argument 
is not valid, since the superfluidity breaks down due to 
the dynamical instability when the effective mass is neg- 
ative. 
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